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The huge data sample accumulated at the KEKB storage ring allows for dedicated analyses in charm spec- 
troscopy. A new, narrow resonance with a mass of 3.872GeV was found in decays of B mesons. Its properties 
remain unexplained. Results on various properties of the D s j(2317) resonances have been updated. Older re- 
sults on the unexpectedly large cross section for double charmonium production in e + e~ annihilation have been 
confirmed and a refined analysis is presented. 



1. Introduction 

Although known as a B-factory, the asymmet- 
ric storage ring KEKB enables the Belle detector 
to explore the charm sector. Charmed hadrons 
are not only the dominant decay product of b 
flavoured mesons, but, in addition, the decay 
kinematics give valuable constraints allowing the 
detection of previously unobserved particles and 
the determination of their quantum numbers. 

Due to the excellent performance of the KEKB 
accelerator, the data sample available for charm 
studies at Belle is unprecedented. Over a year 
ago, KEKB reached its design luminosity of C = 
10 _34 /cm 2 /s and since then shifted its instanta- 
neous peak luminosity to 13.92nb _1 s _1 . Until 
July 2004, Belle accumulated a total integrated 
luminosity of about 287 fb _1 . About 257 ftr 1 
was recorded at the T(4S) resonance, with the 
largest part of the remaining 29 fb^ 1 about 60 
MeV below the resonance. In the following a few 
recent results from the Belle collaboration based 
on smaller data sets are presented, see £Q and ref- 
erences therein for more details on the analyses. 

2. X(3872) 

In decays of B mesons, the b — > ccs transition is 
of great importance due to the CKM matrix ele- 
ments involved. Both vertices of the exchanged 
W boson carry the largest possible matrix ele- 
ments for this decay. So, final states with double 
charm and charmonium states can be produced to 



a vast amount. For charmonium, B decays turn 
out to be a complementary testing field to the 
production of these states in e + e~ annihilation 
at lower energies. 

In general, two variables are commonly used 
for selecting fully reconstructed B candidates, i.e. 
candidates which while decay chain has been com- 
pletely identified and reconstructed. Both make 
use of the fact that at this center-of-mass energy 
(CME), both B mesons are produced back-to- 
back with a well known energy, E — \fs/2. One 
is the beam constrained mass Mb c , where for the 
invariant mass of the B candidate, the energy 
of the candidate is replaced by the well known 
energy E : Mb c = \/E 2 — 0i and - The other one is 
the energy difference AE = E — E can d between E 
and the energy of the B candidate. 

Selecting B ± — > J /^n + TT^ K ± decays by re- 
quiring Mbc of the B candidate to be within 
5.271 GeV < M bc < 5.289 GeV and the energy 
difference AE to be within |AE| < 30 MeV, a 
clear peak at around 600 MeV in the mass dif- 
ference between the J/vP plus the tt + tt~~ system 
and the J/vfr can be identified with the (3770). 
A second peak about 170 MeV higher than the 
one from the \E'(3770) can be seen. MC stud- 
ies showed that this could not be produced by 
e.g. reflections of any known particle, see Fig. 
^ It has therefore been attributed to a new 
particle, dubbed the X(3872), as it's measured 
mass is 3872.0 ± 1.0 MeV 0. It's resolution 
has been determined to be less than the detec- 
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tor resolution for this channel of 2.3 MeV. The 
product branching ratio for J5 ± — > X(3872)_fC ± 
has been determined to B(B ± -> X(3872) J ftT ± ) = 
(1.3±0.3) x 1CP 5 . Its large value suggests that the 
discovery channel is one of the, if not the domi- 
nant decay channel of the X(3872). 



Exp. mass [MeV] width [MeV] 

Belle 2 3872.0 ±0.6 ±0.5 <^3 

CDF 3871.3 ±0.7 ±0.4 < det. res. 

DO @j 3871.8 ±3.1 ±3.0 < det. res. 

BaBar 3873.4 ±1.4 

average 3871.66 ± 0.45 ± 0.31 
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Various two particle invariant masses do not 
exhibit a noticeable structure except the invari- 
ant mass of the tt + tt~ system, see Fig. |2l Similar 
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Figure 1. The mass difference between the 
J/'I , 7r+7r~ system and the J/\& itself. The large 
peak at 600 MeV is due to the decay of the 
#(3770). The second, smaller peak at 770 MeV 
cannot be explained by any known particle and 
has been attributed to a new particle, called the 
X(3872). 



Subsequently, the X(3872) was confirmed by 
CDFg], D0@] and BaBarjSj. Below, the mass 
and width measurements of the X(3872) for all 
four experiments are listed, together with its 
weighted average calculated from the first three 
experiments only as the BaBar result quotes only 
a total uncertainty. The DO measurement has 
been transformed from AM(7r+7r~ J/* - J/*) to 
M(X(3872)) using the current world average of 
M(J/*). 
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Figure 2. The mass distribution of the tt + 7t~ sys- 
tem clusters at large values close to the kinematic 
limit. This suggests, it could originate from a p°, 
which would constrain the C-parity of the new 
particle X(3872). The hatched histogram shows 
data from the sidebands. 



to the corresponding distribution for the 'I' (3770), 
the dipion mass prefers larger values, close to the 
kinematic limit of around 770 MeV. It's distribu- 
tion suggest, the dipion system could originate 
from a p . This would have the direct conse- 
quence of constraining the C-parity of the new 
particle. Would the dipion system actually be 
a p°, the X(3872) must be a C = +1 state as 
C = — 1 particle cannot decay into two C = — 1 
daughters. Furthermore, in this case, the decay 
into X(3872) -> J/*7r°7r° would be forbidden. 
Results on this decay mode can be expected once 
more data is accumulated. 

A second noticeable fact is the absence of the 
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decay of the X(3872) into DD, a dominant decay 
for the considerably lighter 'J (3770). In decays of 
the form B — > DDK the following limit has been 
determined to be 



T(X(3872) -> £>£>) 



< 7 



(1) 



T(X(3872) -» tt+tt-J/*) 
where the corresponding ratio for the '5(3770) is 
r(*(3770) -» £>£>) 



r(*(3770) 



> 160. 



(2) 



This disfavours strongly an assignment of natural 
quantum numbers like ++ or 1" to this new 
state. 

As the X(3872) decays into a charmonium 
state, it is plausible to consider it to be a char- 
monium as well. Grouping all reasonable char- 
monium candidates, i.e. candidates with a sig- 
nificantly large branching ratio into J /^t: + tt^ , 
into groups depending on their C parity leaves 
3 candidates for each possible parity assignment. 
These are the h! c , \&2 and \&3 for the C = —1 as- 
signment and the r)", x' c i an d Vc2 for the C = +1 
assignment. In the following, it will be shown 
that the measured properties of the X(3872) dis- 
favour all six possible charmonium assignments 
listed above. As for the \£ 3 , a large branch- 
ing fraction into DD* is expected. However, 
it has been argued, that the large angular mo- 
mentum strongly suppresses this decay mode. 
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2.1. Is C(X(3872)) = -1 ? 

One obvious candidate for the X(3872) is the 
h' c , for which a sizable branching ratio into 
7r + 7r _ J/\I r is expected and its predicted mass is 
rather close to the mass of the X(3872). Here, 
a distribution worthwhile to check is the angu- 
lar distribution of the in the rest frame of the 
X(3872) with respect to the flight direction of the 



J/^P in the rest frame of the X(3872) candidate. 
For the h' c as well as any other l + ~ state, the dis- 
tribution will be proportional to sin 2 9. However, 
the measured distribution shows clearly a differ- 
ent behaviour as indicated by the reduced x 2 of 
X 2 /d.o.f. = 75/6. 

As for the ^ and ^3, again, the measured 
7r + 7r~ distribution does not agree with the ex- 
pected distribution for these two particles. Addi- 
tionally, for these particles, radiative decays into 
7Xd and 7Xc2 would have branching ratios more 
than two and a half times larger or more than 
three and a half times larger than that of the 
discovery mode, respectively. However, no signal 
has been found and upper limits on the ratio of 
branching ratios have been determined to be 

X(3872) - 7Xci 



< 0.9 



(3) 



< 1.1 



(4) 



X(3872) -> tt+tt-J/* 
and 

X(3872) -> 7 xc2 
X(3872) -» tt+tt-J/* 

respectively. This rules out also these two assign- 
ments. 

In conclusion, no good charmonium candidate 
with C(X(3872)) = -1 is left. 

2.2. Is C(X(3872))=+1 ? 

This would mean that the ir + n~ system is a p°. 
However, the decay would be parity-violating and 
therefore strongly suppressed. Why the discovery 
channel has such a large branching ratio, would 
still have to be explained. 

Here, the candidates are the tj", x'd an d Vc2- 

The mass of the 77" is expected to be in the 
order of <~ 4000 MeV, as it should be about ap- 
proximately 40 MeV below its partner, the %f(3S) 
with m(*(35)) = (4040 ± 2) MeV. 

For the XciJ the radiative decay into jJ/^> is 
expected to dominate over the discovery channel. 
No signal is observed, strongly disfavouring this 
possible assignment. 

For the T] C 2, it is expected that the decay into 
Vc2 — » rj c n + TT~ dominates over the discovery 
channel. However, the large branching ratio of 
the discovery channel would result in a extremely 
large branching ratio for this decay mode. This 
makes this assignment also very unlikely. 
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2.3. New Type of Matter? 

One striking fact however remains: The mass 
of the X(3872) is within the small uncertainties 
identical to the sum of the masses of the D° and 
the D*°. Many years ago such molecule- type 
bound states of mesons were predicted, see e.g. 
references in j2j. 

With more data, the nature of this new particle 
will be more explored. 

3. D sJ (2317) and D sJ (2457) 

Over a year ago, more particles were discov- 
ered, which didn't fit well into the existing quark 
potential models. In c + e~ annihilation, BaBar 
[5] claimed the discovery of a narrow state with 
a mass of 2317 MeV decaying into D+tt° and 
an unexplained excess in D+7r°7 around a mass 
of 2460 MeV. Later, CLEO |7| confirmed the 
first peak and established the second as another 
new particle decaying into D s j(2457) — * D+ 7r°, 
D+ — > D+7. Belle then confirmed both particles 
and identified them also in decays of B mesons 

003 

According to existing quark potential models, 
the masses of these two particles should be more 
than 50 MeV higher than the missing meson 
states containing cs, the 3 Pq and the admixture of 
the 1 P\ and the 3 Pi. The current world averages 
of these two new states are (2317.4 ± 0.9) MeV 
and (2459.3 ± 1.3) MeV, respectively. 

Except for the too low masses, the other prop- 
erties are in a good agreement with the expecta- 
tion for these two cs candidates. 

As both particles are too light to decay into 
DK and D*K, respectively, their dominant de- 
cay, assuming they are the missing cs states, will 
be the modes in which they have been discov- 
ered. Furthermore, the width of radiative and 
other decays are with the current, limited statis- 
tics, well in agreement with the theoretical ex- 
pectations of Bardeen et al. ^Hj an d by God- 
frey Only one upper limit by CLEO on the 

ratio of Z < 0.06 cannot be explained 

by above mentioned models. The only other up- 
per limit on this ratio by Belle determined in de- 
cays of B mesons, is significantly weaker. Ex- 
cept for the discovery channels, only two other 



decay modes have been observed, both for the 
D sJ (2457). These are D sJ (2457) -> D+ 7 seen by 
Belle in both B decays (0.38 ± 0.13) and in e + e~ 
annihilation (0.55 ± 0.15) as well as by BaBar in 
B decays only (0.38 ± 0.12). Latter one is the 
weighted average over various B decay modes. 
The other observed decay mode of the D s j(2457) 
is the decay D s j(2457) — > D+7r + 7r _ with an rel- 
ative branching ratio of (0.14 ± 0.05) w.r.t. to 
discovery channel, observed in e + e~ annihilation 
by Belle. This number slightly, i.e. with a signif- 
icance of about 1 er, disagrees with the two upper 
limit set by Belle in B decays (< 0.10) and by 
CLEO (< 0.08). Other upper limits on decay 
modes are also in the (10-60)% order, and agree 
rather well with the predictions. 

For the lighter D s j(2317), only predictions for 
the decay D s j(2317) — > D+ 7 exist. Here, the 
prediction by Bardeen is closer to the stringent 
limit by CLEO mentioned above, the prediction 
by Godfrey is higher by a factor of three. Ex- 
perimentally, two more upper limits for decays 
of the D s j(2317) have been determined, however, 
no theory predictions for these two modes exist. 
These two decay modes are D s j(2317) — > D+7 
and D s j(2317) — > D+7r + 7r _ , which are forbidden 
for the missing cs states. Stringent limits in the 
percent or even the permille level support the hy- 
pothesis, this particle is in fact the so far unob- 
served cs state. 

More information about the spin of an parti- 
cle can be derived from angular distributions like 
hclicity distributions, decay angle distributions in 
the rest frame of the mother particle w.r.t. the 
flight direction of the mother in the rest frame of 
the B. 

For the D s j(2317) the helicity angle distribu- 
tion for the decay D s j(2317) — > D+7T is in excel- 
lent agreement with the spin hypothesis . 

Also for the D s j(2457), the helicity angle dis- 
tribution in D s j(2457) — > D+ 7 supports the spin 
1 assignment of this particle, likely being the so 
far missing cs state . 

In conclusion, the new discovered particle, 
named the D sJ (2317) and the D sJ (2457), be- 
have like the unobserved cs states, however their 
masses are significantly too low. 
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4. Double cc production 

Another very nice and elegant analysis show- 
ing surprising results is the investigation about 
the double cc production. Here, a J/\& is recon- 
structed in the very clean lepton channel, then, 
simply, its recoil mass is plotted. The effects of 
ISR photons are taken into account, see |12I13| for 
details. A momentum cut removes contribution 
from decays of B mesons. 





L jM 


\ J, ^jj- 


Wr 





2.2 2.6 3 3.4 3.8 

M re coii(JM ( QeV/c2 ) 



Figure 3. The recoil mass spectrum against 
a fully reconstructed J/'J' in e + e~ annihilation 
events. A clear enhancement of activity above 
the cc threshold is visible, in contrast to predic- 
tions by NRQCD. 



The recoil mass spectrum does not reveal a 
large signal until about 3 GeV, the charm pro- 
duction threshold, where the signal strongly in- 
creases. 

Zooming into the threshold region, see Fig. 
reveals three clear distinct peaks, in a previous 
analysis |12| identified with the rj c , XcO and the 
7] C (2S). The masses of these particles were con- 
sistent with the world average. For the ?7 C (2S) , 
it was only the second observation. Note, that 
e.g. the production of two J/W from a single 7 is 
forbidden by conservation of C-parity. 



Another surprising result of the old analyses 
were that there is almost no activity below 3 GeV, 
but a very strong activity above. 

So far, the only theory predicting cross sec- 
tions of charmonium production at this center- 
of-mass energy is non-relativistic QCD, NRQCD 
for short. It predicts the 3/^gg production to 
be dominant with a total cross section of about 
~ 1 pb. The colour singlet production of J/^g 
might dominate the phase space as a quasi two 
body decay at the endpoint with up to ~ 0.5 pb. 
Double cc production via gluon splitting into cc 
is expected to be suppressed by about a factor of 
10-20 w.r.t. to J/%3. 

The measurement however finds the J/^gg 
process to be negligible and the double cc pro- 
duction dominates with about 80%. The total 
cross section for double cc production is about 
one order of magnitude too large compared to the 
prediction of NRQCD. 

Several, alternative explanations arose, e.g. 
that under the ?y c peak, significant contribution 
from c + c~ — > 77 — * J/^J/ty could be hidden. 
Glueballs could interfere with the charmonium 
states and increase the cross section, or, a bias 
in the momentum scale can deteriorate any inter- 
pretation, that the first peak would no be due to 
the r) c . 

In a refined analysis [?] performed on a larger 
data set of 155 fb _1 , the fit to the recoil mas 
spectrum now allowed for all narrow resonances 
in this mass region. To the r) c , XcO, ?7c(2S) allowed 
in the previous analysis, the J/^, Xci(2) and the 
^(25) were added. However, the fitted yield for 
the last three particles was consistent with zero, 
even negative for the 3/^> and the \l/(25). As 
no signal was observed, the mass and width of 
these particles were fixed to their world average 
values. For the other three particles, these values 
were floated in the fit. The fitted yields are listed 



below 


cc 


mass [MeV] 


N 


signif. 


Vc 


2972 ± 7 


235 ± 26 


10.7 


J/* 


fixed 


-14 ±20 




XcO 


3407 ± 11 


89 ±24 


3.8 


Xcl(2) 


fixed 


10 ±27 




%(2S) 


3630 ± 8 


164 ±23 


5.3 


*(2S) 


fixed 


-26 ±29 
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In order to confirm the peak around 3 GeV is 
caused by ry c and not by J/ 'J, various tests have 
been performed. First, about 230 signal events 
in the ry c mass region allow fully reconstruct a 
few events. The ry c was reconstructed in the rj c — > 
KgK + 7r~ and the r\ c — > 2(K + K~) decays. From 
Monte-Carlo, 2.6 events were expected to be re- 
constructed. In the data, 3 events have been ob- 
served, in excellent agreement with the expecta- 
tion. Additionly, in the dilepton channel of the 
J/^> no event events has been observed, in accor- 
dance with the expectation. 

Furthermore, an angular analysis of the pro- 
duction and helicity angle distribution for the 
reconstructed J/ "J have been performed for all 
events in the three signal peaks. For the produc- 
tion via e + e~ — > 7* — > J/^Pcc, due to conser- 
vation of angular momentum of the exchanged 
7*, the production and the helicity angle dis- 
tributions have to be equal. For processes like 
e + e~ — ► 77 — ► J/^cc no such constrain exist. 

All distributions have been fitted with a func- 
tion of the form f(x) ~ 1+acos 2 with a floated. 
Within large uncertainties, the results are consis- 
tent with the hypothesis of a 1 7 exchange. The 
last column in the table below actually assumes 
the production via the 1 7 exchange, reducing 
slightly the uncertainty on a. Also listed in the 
last row arc the predictions for a glueball. 
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Additionly, the cross sections for the produc- 
tion of double charmonium have been determined. 
The results are listed in the table below. Note, 
that for the decay of the charmonium accompany- 
ing the J/^, two or more charged particles were 
required in the analysis, reducing the visible cross 
section and therefore enhancing even more the 
discrepancy between theory prediction and exper- 
iment. 



cc cTBom x B(cc^> 2charged) [fb] 



Vc 


25.6 ±2.8 ±3.6 


J/* 


< 9.1 


XcO 


6.4 ± 1.7 ± 1.0 


Xcl(2) 


< 5.3 


r/c(2S) 


16.5 ±3.0 ±2.5 


9(2S) 


< 13.3 



5. Conclusion 

Over the last years, the 6-factories made sev- 
eral new unexpected discoveries, mostly involv- 
ing charmed mesons. An apparent charmonium 
state, the X(3872), does not fit at all into the 
charmonium spectrum. It might represent a new 
form of matter, a molecule build out of mesons. 
Other new discovered particles in the cs system 
behave like expected, however, their masses are 
significantly too low. And last, but not least, the 
production cross section of double charmonium in 
c + c~ annihilation is still about an order of mag- 
nitude larger than predicted by theory 
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